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Abstract 

This paper discusses the problem of gravitational perturbations of ra- 
diating spacetimes. We lay out the theoretical framework for describing 
the interaction of external gravitational fields with a radiating spacetime. 
This is done by deriving the field perturbation equations for a radiating 
metric. The equations are then specialized to a Vaidya spacetime. For 
the Hiscock ansatz of a linear mass model of a radiating blackhole the 
equations are found separable. Further, the resulting ordinary differential 
equations are found to admit analytic solutions. We obtain the solutions 
and discuss their characteristics. 



1 INTRODUCTION. 

The study of gravitational perturbations can be traced back to the famous 
Einstein-Infeld-Hoffman paper of 1938 M which pioneered the treatment of the 
two body problem in general relativity. In 1957 Regge and Wheeler Q addressed 
the problem of the stability of a Schwarzschild black hole. Later, in his study 
of perturbations of a rotating black hole ( |!| and later papers), Teukolsky was 
able to put the discipline on a stronger footing. However, little progress has 
been made at extending this success to cover the radiating cases. The problem 
of perturbing a radiating spacetime with integral spin fields has not received 
the attention it deserves. This, despite the fact that most astrophysical objects 
radiate. From regular stars to supernovae, from quasars to primordial black 
holes one finds that the inhabitants of our universe are generally non-static. 

In the present paper, we develop a framework for discussing the problem of 
how external gravitational fields may interact with radiating spacetimes. This 
is done by deriving the field perturbation equations. It is found that two such 
equations are sufficient to describe all the non-trivial features of the perturbing 
gravitational field. We find that one of these equations decouples completely 



and is homogeneous in one of the field components. The equation for the other 
field component contains, in its source terms, several perturbed and therefore 
undetermined quantities. Using a systematic approach we are able to determine 
all these quantities completely, in terms of the former field component. The re- 
sult is that all the perturbations are described by only two field components 
which satisfy two partial differential equations. The equations are then special- 
ized to a Vaidya spacetime. For a particular model of a radiating black hole the 
equations are found separable. Interestingly, the resulting ordinary differential 
equations are found to admit analytic solutions. We obtain these solutions and 
discuss their characteristics. 

The mathematical framework used in this paper is the null tetrad formalism 
of Newman and Penrose (hereafter NP formalism) j(| . In section II we give a 
brief description of the background geometry, the radiating spacetime of Vaidya. 
In section III we derive the perturbation field equations for a general non- vacuum 
type D spacetime and adapt these equations to the Vaidya spacetime. In section 
IV we calculate the perturbed quantities in the source terms to arrive at the final 
working field equations. It is demonstrated in section V that these equations are 
separable for the Hiscock linear model |0| of a radiating black hole. In section 
VI we obtain and discuss some of the solutions and we conclude the discussion 
in section VII.. 



2 THE VAIDYA SPACETIME. 

2.1 The Metric 

In this analysis we perturb the Vaidya spacetime j|] with incoming external 
gravitational fields. The Vaidya geometry, the simplest of the radiating space- 
times, is non-rotating and spherically symmetric. The energy-momentum tensor 



T, 



pk^k v 



describes a null fluid, (k^k^ = 0), of density p with radial flow, k 2 = k 3 = 0. 
Using this energy-momentum tensor to solve the Einstein field equations one 
obtains |Io| a line element, in retarded coordinates, given by 



ds 2 = 



1 - 



2m (u) 



du 2 + 2dudr~r 2 dfl 2 . 



(1) 



Here, u is the retarded time coordinate and m(u), the mass, is a monotoni- 
cally decreasing function of u. 

It is convenient to introduce a null tetrad basis z mfi = (7 M , n^, m^, fhp) at 
every point in this spacetime. The metric tensor then becomes 
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where r\ mn is the flat spacetime metric. Following Carmeli and Kaye 
choose the covariant form of the null tetrad basis as 



we 



h = 5 n 



1 - 



2m (u) 



-—[S^+ism 66*] 



(3) 



sin 65*] 



The contravariant vectors, z m , considered as tangent vectors, define the direc- 
tional derivatives as 



y — 7 



(4) 



In the Vaidya spacetime these directional derivatives are given (from equa- 
tion p|) by 



D 
A 

5 
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1 - 



2m (u) 
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(5) 



One finds that the only surviving spin (Ricci rotation) coefficients (8| of the 
Vaidya spacetime are 



P 



2V2 ; 



cot 9, f3 = —a, 



1 

27 



2m (u) 



and 7 



m (it) 



(G) 



The only surviving tetrad component of Ricci tensor is 



3 



-m (u) 



(7) 



where m = f 1 . With this component, it is easily shown that 



Tfj, v — 2$22^/j^ 



771 (u) 



(8) 



Further, the only non-vanishing tetrad component of the of the Weyl tensor 



is 



*2 = 



-771 (u) 



(9) 



The Vaidya spacetime is then || said to be Petrov type D with repeated 
principal null vectors Z M and n M . The three optical scalars are found, to be 
ct = cj = 0, 9 = —-. The metric contains two shear- free, twistless and 
diverging geodetic null congruencies. 



3 The Perturbed Field Equations 

3.1 The Type D Spacetime. 

In this section we develop the gravitational field perturbation equations for a 
general Petrov type D spacetime. We start with the Newman-Penrose (NP) 
equations. The full set of the NP equations can be found in several publica- 
tions, for example |^|. We only mention here that the set is made up of first 
order differential equations which, in the NP formalism, replace the Einstein 
field equations. The equations link together the null tetrad basis, the spin 
(Ricci rotation) coefficients, the Weyl tensor, the Ricci tensor and the scalar 
curvature. In using this formalism to do perturbation analysis one first spec- 
ifies the perturbations of the geometry. Here we shall write the tetrad of the 
perturbed spacetime as 



l = l b + F, n = n b + n p , m = m b +m p , m = m b +m p . (10) 

where the superscripts b and p refer to the background and the perturbing 
quantities, respectively. Since all the other field quantities are expressible in 
terms of the tetrad m their perturbed forms can be written down. For example: 

We shall, in general, assume that the perturbations in the basis vectors are 
sufficiently small so that only their first order contributions are significant. The 
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field equations are, then, first order in the perturbing fields, linearized about 
the background quantities. 

In type D spacetimes one can choose the tetrad vectors (see section II) 1 and 
n so that certain spin coefficients and certain Weyl scar components 



5*3 - D^ 4 + <5$2i - A$2o = 3A* 2 - 2 (a + 2tt) * 3 + (4e - p) * 4 

-2i/$io + 2A$ U + (27 - 27 + p) $20 + 2 (f - a) $21 - ^22, (12) 



A* 3 - (5*4 + <5$22 - A$2i = 3i/* 2 - 2 (7 + 2/i) * 3 + (4/3 - r) * 



-2i/$u - i>$20 + 2A*i2 + 2 (7 + p) $21 + (f - 2/3 - 2a) $22, (13) 
from the spin coefficient system of equations we take 

AA - 5u = - (fi + p) A - (37 - 7) A + (3a + /3 + tt - r) v - $4. (14) 

We have complete knowledge of the geometry of the background spacetime, so 
we write down the equations only in those terms that make first order con- 
tributions to the perturbed field quantities. Making use of equations [ll] on 
equations |l2|, [l3] and [l4|, respectively, it is seen that for a perturbed Petrov type 
D spacetime 




and 



3A P * 2 + {S + 2a + 4tt) 



(D + 4e - p) 




(15) 



-3i/ p * 2 + (A + 2 7 + Ap) * 3 P - (4/3 - t + 8) $r£ 



(2 7 + 2/2 + A) 
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+ (-8 + f - 2/3 - 2a) p $ 22 + (-5 + f - 2/3 - 2a) , (16) 

and 

(A + 3 7 - 7 + fi + p) X p - (6 + 3a + + tt - f ) v v + * 4 P = 0. (17) 

We note that the perturbed equations are coupled both in the Weyl tensor 
components and the Ricci tensor components. They also contain unknown spin 
coefficients and directional derivatives. In attempting to decouple the equations 
above we use an approach akin to Teukolsky's H. After some algebra we find, 
on eliminating $3 between the equations [l^ and |l6|, we are left with 

[(A + 37 - 7 + An + p) {D + 4e - p)] 

- [(5 + 3a + fi + 4tt - f) (4/3 - r + S) - 35- 2 ] = Q 4) (18) 
Now, under the interchange 1 ^ n, m ^ m, the full set of the NP equations 



is invariant 1 1 . This symmetry is not destroyed by the choice of 1 and n which 
gave equations O. One finds under this interchange that JtJ 



k <=* — v, p <=^ — p, a —A, a — /3, e <=^ — 7 and 7r <=± — f , (19) 
Applying this to equation [l^ we obtain the following equation for 



(D - 3e - e - 4p - p) {6 - 3/3 - a - 4r + tt ) (5 - 4a + tt) 



- (D - 3e - e - 4p - p) [(A - 4 7 + /x) - 3* 2 ] *<? - Qo- (20) 
Here, the source term Qo is given by 

Qo = (6 - 3/3 - a - 4r + tt ) (D - 2p - 2e) $ | 
- (£ - 3/3 - a - 4r + tt) (£ + tt - 2a - 2/3) 
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+ (8 - 3/3 - a - 4r + tT) (S + tt - 2a - 2f3) p $ 



(id 



(D - 3e + e - 4p - jo) (<5 - 2/3 + 2tt) $, 



- (£> - 2e + 2e - p) + (£> - 3e + e - 4p - p) a v (A - p) $ 00 . (21) 

Equations [l8] and |2^ describe the gravitational field perturbation equations 
for a general type D spacetime with sources. 

3.2 The Perturbed Vaidya Spacetime. 

We shall now adapt the equations above to the problem of perturbing a Vaidya 
spacetime with gravitational fields. The important features to deal with here 
are the perturbations in the source terms Qo and Q4 of equations ^| and [2^. 

The energy-momentum tensor associated with the radiation in the Vaidya 
spacetime is usually interpreted via geometrical optics. Carmeli and Kaye Q 
have, in fact, shown that the associated radiation field which has a monopole 
structure can not be identified as a source-free electromagnetic field. There- 
fore in addressing the problem of perturbing this radiation, one can only dis- 
cuss direct perturbations on the energy-momentum tensor. Now, this energy- 
momentum tensor is quadratic in the metric p7| . It follows from these con- 
siderations, then, that the lowest order perturbations in the scalar components 
&mn will be quadratic in the metric perturbations. However, we assumed from 
the beginning that the perturbations in the basis vectors are sufficiently small 
so that only their first order contributions are significant. In our linear theory 
we shall, therefore, disregard these perturbations in the scalar components <fr m n 
contributions. Nevertheless, we still have various perturbed spin coefficients and 
directional derivatives to deal with in the source terms. 

It is recalled (see section II) that the only quantities which survive in the 
Vaidya spacetime background are: the spin coefficients a, j3, 7, p, /i, the Weyl 
tensor component ^2 and the Ricci tensor component $22 • Moreover, all these 
spin coefficients are real and a — —j3. On applying the observations made above 
to the field perturbation equations Il8| and |2C| we find that 



[(A + 2 7 + 5/i) (D — p) — (5 — 2/3) (8 + 4/3) - 3* 2 ] = Q' 4 , (22) 
[(£> - hp) (A - 4 7 + p) - (8 - 2/3) (6 + 4/3) - 3* 2 ] * P = 0, (23) 
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where the source terms are now given by 

Q' 4 = [(A + 2 7 + 5/x) a p - (6 - 2(3) (S p - f p + 2~(3 P + 2a p ) + X p (D + p)} $ 22 , 

(24) 

and clearly 

Q'o = 0. (25) 

Note that in equation ^4] the derivative operators A and D must act on $ 22 . 

We now have two equations, |22| and ^3| which describe the two Weyl tensor 
components \1/ P and 'F/. It can be shown that these two components are 
sufficient to describe all the non-trivial features of the perturbing fields. The 
proof of this sufficiency is achieved by showing || that only ^ p and 'F/ are 
invariant under gauge transformations and infinitesimal tetrad rotations. This 
invariance, in turn, makes them completely measurable physical quantities. 

The source terms Q' 4 in equation |2^ still contain (see equation ^4|) pertur- 
bations in as many as seven different quantities. In the next section we shall 
derive equations that describe these unknowns. In doing so, we shall discover 
that those perturbed quantities that do not vanish are all completely describable 
in terms of the one of the two field components. 



4 The perturbed quantities in the source terms. 

The aim of this section is to calculate, (in terms of the known background quanti- 
ties and the fields *k p and ^ p ) the perturbed quantities, <& p 2 , ct p , (3 P , i~ p , <r p , X p 
and 5 P which appear in the source terms Q\ of equation ^2], (see also equation 
p4| ). It is worthwhile to point out that the system of equations we have con- 
structed in the previous section must be consistent with the freedom we have 
in the choice of both the tetrad frame and the coordinates. In particular |Q, 
from the 6-parameter group of homogeneous Lorentz transformations, we have 
six degrees of freedom to make infinitesimal rotations of the local tetrad-frame. 
Further, we have four degrees of freedom to make infinitesimal coordinate trans- 
formations. Thus, altogether we have a total of ten degrees of gauge freedom. 
We are free to exercise these available degrees of freedom as convenience and 
occasion may dictate. 

As has been shown ||, in a linear perturbation theory, and \I/ 4 P are 
gauge invariant while ^ p , ^ p and , F 3 P are not. This means that we can subject 
the tetrad null-basis to an infinitesimal rotation in which ^ p and <F 3 P vanish 
without affecting ^/ p and \& p . In making this choice of gauge, we have used 
up only four, out of ten, degrees of freedom. It can be shown (see for example 
[fr] Chapter 1, equations 342 and 346) that in linearized type 1 rotations both 
\E , 1 P and $ p 2 vanish while in type 2 rotations both do not. It follows, then, 
that a gauge in which both ^> p , 5 , 3 P and $ p 2 vanish can be chosen. Further, in 
this gauge, ^ p clearly vanishes. Under these circumstances, the NP equations 
show that the linearized Bianchi identities take on a simpler form. For our 
purposes, the equations we need to consider from this set are: 
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A* p = (4 7 - fj,) * p + 3o** 2j (26) 

-3(5 p * 2 = -9t p * 2 -2 k p $22, (27) 

-£>tf 4 p = 3A p tf 2 - ptf/ - a p $ 22 , (28) 

= -k p $ 22 , (29) 
4.1 Calculation of ct p and A p . 

From equation ^6] and the fact that all the background quantities here are real 
we find that 



+ (30) 



3* 2 

Further, using equation au on equation GH gives 



We note that the perturbed spin coefficients A p and a v display a definite 
dependence on vE^ and/or \1/ P . In the rest of this section we shall derive ex- 
pressions for the remaining perturbed quantities. 

4.2 The perturbation matrix for the basis vectors. 

In order to determine the perturbed quantities t p , a p , (3 P and 5 P and their 
relations , it is necessary to study the effects of the perturbations on the basis 
vectors, (Z M , n M , m/ i ) m* 4 ). For compactness, it is convenient to introduce the 
following index notation: 



i 1 = I**, Z 2 = n**, i 3 = m**, I 4 = m**. (32) 

We can write ]7j the perturbations (i — 1,2,3,4) in the vectors, as 

linear combinations of the unperturbed basis vectors Z*. Thus 



l (p)i = pi l 3 i (33) 
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pi 


p 1 

r 2 


P 1 
r 3 


P 1 


p\ 


p2 
r 2 


P 2 
r 3 


P\ 


p'i 


p3 


p3 


P\ 


pf 


pi 
r 2 


P 4 


pf 



where the P* are elements of a matrix P that describes, completely, the pertur- 
bations in the basis vectors. Explicitly, 



(34) 



The I 1 and I 2 are real while the / 3 and I 4 are complex conjugates. It follows, 
then, that the matrix elements P\, P\, P\ and P\ are real while the remaining 
elements of P are complex. Moreover, the elements in which the indices 3 and 
4 replace one another, are complex conjugates. For example, P\ = (Pf) 

4.3 Perturbations in the Angular functions, 5 P , t p , a p , and 

(3 P . 

The perturbations in the directional derivative S p , are given from equations^ 
and t33L by 



F = (PV 4 ) P . 

But from equations and p4 we see that 



(35) 



=PU j = P\l X + Pi? + Pjl 3 + P\l 4 



Using equation 36 on Equation [h3j shows that 



(36) 



PjD + P 4 2 A + P 4 3 S + P 4 ~6. 



Thus, if we operate with S p on the background \&2 we get 



(37) 



F<f 2 = P\D^ 2 + P\^ 2 + Pt^2 + PiSV 2 - 

Now, in the background 



(38) 



*2 = 



— m{u) 



Substituting for ^2 in equation 38 and using definitions of the operators in 
equations |^ we find that 



F* a = P f\^pi + Pl hil\ -P 4 Jl- ^1 . (39) 
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Recall that in the Vaidya spacetime background, the component $22 = 
— m ^ is non- vanishing. Using this on equation ^9] shows that n p must van- 
ish. It follows, then that equation |27| becomes, 



<f $ 2 = 3f p ^ 2 , (40) 

We see immediately, that the results expressed in equation ^ will be incon- 
sistent with the eigenvalue equation ^ unless the P\ vanish, so that 

J3\T, _ D 4 [3™ («)] _ r>4 ( 



= ^T^F 1 = -3Pf " *2- (41) 



Equations [lO] and then, show that 



-Pi(-). (42) 



Moreover, equation ^ along with the condition that the P\ vanish means that 
whenever S p acts on a function with no angular dependence its only contribution 
is 



h> ,,a d _„ d 



5 =^0-r = - fPr 3-r (43) 

We now have, in equation ^3], a general relationship between S p and f p . 

Next, we need to deal with a p and Jf . It is known || that if the null vectors 
1^ are tangent to the geodesies and equal to a gradient field, then 



p = p and t = a + (3. (44) 

These conditions are fulfilled in all Type D space-times. In particular, the 
unperturbed Vaidya space-time satisfies 



a + (3 = t = 0. 

Consequently, in our linear perturbation analysis we should have 

( a + P) P = fP. (45) 
Consider, now, the second of the source terms in equation E4] which reads as 



(5 - 2/3) (S p -f p + 2f + 2a p ) $ 



22- 
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Our results in equations ^3] and ^5] when applied to the above expression show 
that 



(S - 2(3) (S p -fP + 2f + 2a p ) $ 22 = 3 (S - 2/3) ? p <f> 2 



(46) 



This achieves the purpose of expressing the effects of the four perturbed quan- 
tities t p , a p , j3 p and 8 P in terms of one of these quantities t p . The next task is 
to express this perturbed quantity in terms of the fields. To this end we utilize 
one of the equations from the spin coefficient set fel . Consider the equation 



St — Act = [per + Xp) + (r + (3 — a) t — (37 — 7) ct — nv + $0 



(47) 



Specializing equation ^ to the Vaidya spacetime and applying our linear pertur- 
bations approach we obtain, on rearranging terms, an equation whose complex 
conjugate is 



(S - 20) f p = (A + p - 2 7 ) g p + P X P . 



(48) 



Equation ^ when substituted into equation [l6] gives (recall S does not op- 
erate on $22) 



(S - 2(3) (S p -f p + 2(3 P + 2a p ) $ 22 = 3 [(A + M - 2 7 ) a p + p\ p ] $ 2 



(49) 



4.4 The final field perturbation equations. 

We are now in position to apply the results of our discussions in this section to 
the source terms of equation ||| as given by equation |||. Using equation |9| on 
equation 24 and rewriting the sources Q\ as Q' we obtain, 



[(A + 2 7 + hp) a p - 3 (A + p - 2 7 ) a p - ip\ p + \ p (D + p)] $ 2 



(50) 



We see that in the above equation the terms in X p cancel yielding the simpler 
result, 



Q' = -2 (A - 47 — p) ct p $ 22 . 
Substituting for a p from equation pfj we find that 



(51) 



Q' = -2 (A - 4 7 - p) 



22 



A - 47 + — ; 



3* 2 



(52) 
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Equation p32| forms the result of our analysis in this section. All the perturbations 
in the sources have now been expressed in terms of the perturbed field only. 
The working field equations (see |2^ and ^) have become 

[(A + 2 7 + 5/*) (D-p)-(5- 2/3) (S + 4/3) - 3* 2 ] = Q', (53) 

and 

[(D - bp) (A - 4 7 + m) - (5 - 2/3) (6 + 4/3) - 3* 2 ] * P = 0, (54) 

where, now, Q' is given by equation [5^. 

Equations ^3], |54| along with ^2] form the main result of our perturbation anal- 
ysis. These equations give the essential features of a gravitationally perturbed 
Vaidya space time. All the non-trivial perturbations are sufficiently described 
by two tetrad scalar components of the Weyl tensor, and ^ 4 P , which com- 
ponents represent the extreme helicity states of the gravitational field. 

We now can rewrite the equations in a form that reveals the dependence of 
the fields on the physical variables of spacetime. Thus, using equations @, [?] 
and ^ on [54|, |5^, and [52| respectively, we find that 



if 



5 9 If 2m (u) \ d 2 3 
dr du r du 2 \ r J dr 2 \ 



m (u) \ d 



d r 



i 

2>2 
and 



2 d 

J^ + ^t9—~2(csc 2 9 + cot 2 , 



d d 
esc 2 9-—, - + 4i esc 9 cot 9—- 

oip z dip 



* * = 



(55) 



d 2 



dudr r du 2 



I i) I / i 2m (u) \ d 2 1 / 7m (u) \ d 



dr 2 



2 ! ^ 4m (u) ' 



1 

Jr 2 



d 2 d 

cot 9 — - 2 (esc 2 9 + cot 2 1 



dr 



89 



d 2 d ' 

esc 2 9 — — x — 4z esc 9 cot 9— — 



where now 



*i = Q 



(56) 



Y = -2 



d 1 / 2m (u) 
du 2 V r 



d 3m (u) 
dr r 2 



1 

27 



m (it) 
v 3m (u) 



d 1 / m (it) 
(9u 2 V 



9r 



m (w) 



1 

27 



(57) 
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5 Separation of Variables 



In this section we seek to separate the equations that were derived in the previous 
section. This separation of variable is effected in two phases. In Phase I we deal 
with the angular variables while in Phase II we deal with the retarded time and 
radial variables. We shall, for now, concentrate on the homogeneous parts of the 
equations. The contribution due to the source terms can always be constructed 
later once a solution for has been obtained. Incidentally, one notices (see 
equation [57l) that the luminosity-mass ratio 3n ^ u ^ , (L — —to) which scales the 
source term will almost always be vanishingly small since for most radiating 
objects the mass being radiated at any given time is much smaller than the rest 
of the body mass. 

5.1 Phase I: The Spin- weighted Angular functions. 

We suppose that the gravitational fields entering the spherically symmetric back- 
ground spacetime are plane waves so that the problem has azimuthal symmetry. 
With this, we then assume that the field equations are separable in the angular 
variables admitting solutions of the form 



*i=(o,4) C",r,0,99) = i= ( O ,4) C"> r >#) e " nv = R P =(±2) (u,r) <Sp=(±2) (0) e. (58) 

Here the subscript p is used to identify a particular spin-s field component by the 
spin weight. For our purposes, the spin weight p only takes on the extreme values 
of ± s corresponding to the extreme helicity states of the field. Explicitly, 
has a spin weight of 2 while $4 has a spin weight of —2. Note, to avoid confusion 
in notation, here and henceforth we discard the superscript p previously used 
to identify the perturbed quantities. 

Substituting [58] in the field equations [55] and |5^ yields the following general 
equation in the angular variables: 



1 d ( ■ n d \ / n , n 2mpcos9 m 2 

-— -— sin 6 — ] + [p - p 2 cot 6 ^5 5 K \S V (9) e mip = 0. 

sin 9 d9 \ d8 V F sin 2 8 sin 2 9 1 



(59) 

This equation along with boundary conditions of regularity at 9 = and 9 = 
7r constitute a Sturm-Liouville eigenvalue problem for the separation constant 
K = p K™. For fixed p and m values, the eigenvalues can be labelled by I. The 
smallest eigenvalue has I = max(p, |m|). For each p and m the eigenfunctions 
pS 1 ™ (9) are complete and orthogonal on the interval < 9 < n, as required by 
the Sturm Liouville theory. In our case, where the background is non-rotating 
the eigenfunctions are, the well known p3|, spin- weighted spherical harmonics: 



?Yr (9, if) = p Sr (0) e mv (60) 
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and the separation constant K is found to be given by 



K= p K, = (l-p)(l+ P +l) (61) 
5.2 Phase II: The radial-null equations 

The separation of variables effected in the last section leaves us with two equa- 
tions for the functions R+ 2 and R- 2 . These functions are coefficients of 2 Y™ (0, <p) 
and -iY™ {6,<p)i respectively, in the spin-2 fields Wo and ^> 2 and are each de- 
pendent on u and r only. On substituting equation ^0| into equations ^5] to ^7] 
one finds that R+ 2 satisfies 



d 2 R 2 (u, r) 5 dR 2 (u, r) 1 / 2m ju) \ d 2 R 2 
dr du r du 2 \ r I dr 



m(u)\ dR 2 (u,r) ( 2 K t - 4) 



d 1 



2r 2 



R 2 (u,r) 



= 0. 



(62) 



and i?_ 2 satisfies 



~ d 2 R- 2 1 dR- 2 1 
dudr r du 2 



2m(u) \ d 2 R- 2 
dr 2 



3 7m(u)\ dR- 2 
r r 2 I dr 2 



+ (^(-^- 4 ) + ^)^(v)=0 (63) 

We deal with equations ||| and |63| separately. First, we shall seek to separate 
equation |6^ in R 2 (u, r). By adopting a change of variables we show that the 
equation is separable for a specific choice of mass function. Thereafter, we shall 
apply this approach to equation |6^. 

5.2.1 Change of variables. 

Equation [32|, as it stands, is not separable. We shall, therefore, find it convenient 
to introduce the following change of variables: let us set 

r 

Then it is seen that 

d_ 

du 



- and £ = ^. 

u r 



(64) 



d 



dr m (t) <9£ 
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d 



dr 2m (t) ' 



(65) 



and 



dr 2 2 [m (r)] 2 d£ 4 [m (r)] 2 d£ 
where, now, m is a function of r but m still means ^° L . 



2 ' 



5.3 Equation for R+2- 



On substituting equations 64 and pq into p2l and rearranging we find that 



-5* 



(e 4 -4e 3 )^+e^ 2 i?2(r,0 = (66) 



where we have set 



ei = (l-p)(l+ P +l) 



(1-2) (1 + 3) 



(67) 



We would like, as an example, to apply our analysis on an evaporating 

is not separable for an arbitrary mass 



66 



blackholc. In general, the equation 
function m (u) . However it can be shown to separate for one particular model 
of such a radiating blackholc. 



Vaidya model for a linearly radiating blackhole. 

In the Vaidya model of a radiating blackhole ]l2]| , the spacetime is, initially, 
Minkowski flat for u < 0. Then at u = an imploding 8— function-like null 
fluid with a total positive mass M forms a blackhole. Hereafter, < u < uq 
negative energy null fluid then falls into the blackhole evaporating the latter 
in the process. One known consequence jjjj is that the spacetime violates the 
weak energy condition. Eventually the blackhole vanishes so that for u > u$ 
the spacetime becomes Minkowski flat again. 

One of the popular models of radiating black holes is the so-called self-similar 
model originally developed by Hiscock, p|. Popular, because from it one can 
construct the quantum energy stress tensor for the entire spacetime. The model 
has been extensively used lately, (see for example |Q and Jlg^. In this model 
the mass is a linear function of the retarded time coordinate u. 
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We shall show, presently, that for the Hiscock linear mass function ansatz 
the above equation is separable. 
Suppose 

( 0, u < ] 

to (u) = I M (1 - A u) , i > u > \ (68) 
I 0, u > i J 

so that 

to = -AM (69) 

where too is the initial mass at u — and A is some positive parameter < A < — 
that scales the radiation rate. 

We shall find it convenient to institute a change of variable u = v + v where 
vq is some fixed value of u. 

too (1 + Xv) , — vq < v < 

m (v) = l m \ " = . 

» to [1 - Aw] , j > v > » 



0, v > 



This seemingly trivial change is important for the following reasons. In our 
problem we would like to discuss the behavior of the gravitational fields in a 
radiating blackhole background. However as we noted above just before u = 
there is no black hole and yet we need the ingoing fields to have been moving in a 
non-minkowski background. This change, therefore, makes it physically possible 
for us to introduce the external gravitational fields into a spacetime that already 
contains the black hole. Mathematically the change makes it possible, as we find 
out soon, to construct complete solutions that include a description of ingoing 
fields. 

Substituting equation p3 into pq and making use of equation p3 we find that 



-5£— h4 



Or - d ^ T )^m(r )+ ,(e-e + Um e) ^ 



- (C 4 - 4C 3 + 16A m e 2 ) ^ + e t e R 2 (r, = (71) 
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5.3.1 Separability. 

We now ask whether the equation |7l] is separable (in the variables r and £) 
admitting a solution of the form 

R 2 (t,0= X 2 (OY 2 (t) (72) 

where X2 and Y% are each functions of one variable, only. Substituting for R2 in 
equation (n] yields two ordinary differential equations: a first order differential 
equation for Y 2 (r), 

wr2 (yM if) = ~ a (73) 

and a second order ordinary differential equation for X2 (£), 



j2 v- j y 

e (e-ti + 4Am ) - € [e 2 - 4C + 4a + 16Am ] + (e, £ + 20a) X 2 (0 = 0. 

(74) 

Here, a is an arbitrary separation constant. Its characteristics are discussed in 
the next section. 



5.4 The function i?_ 2 

Following the same approach as above we find that the equation for the function 
R-2 separates into two ordinary differential equations: 



and 



2 1 dY_2 

m t r 

1 ' Y-7. dr 



d 2 X 



= ~7> 



(75) 



C (r - £ + 4m A) - £ [5f - 4£ + 47] 



+ (e 2 + 2e ; e + 4 7 ) X_ 2 = 0, (76) 

The general characteristics of the separation constant 7 are not different from 
those of a and discussed in the next section. 



6 Solutions 

We have achieved the separation of the original equations ^ and ^ respectively, 
into |7^, [74| and |75| and |7^ for the particular case of a linear mass function. In 
the following sections we shall seek to solve these resulting ordinary differential 
equations and to discuss the solutions in an attempt to draw some physical 
information from them. 
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6.1 The functions Y p (r) 

The first order differential equations for Y 2 (r) and Y_ 2 (r) above can be inte- 
grated immediately. Thus from equation |73| we find that 

Y 2 (v) = exp(-a x J - e ■ ( 77 ) 

V Jo m (l-\v)J 

Now < A < 1 and in fact for most radiating bodies Ac 1. This allows us to 
expand the logarithmic expression In (1 — Ait) in the solution so that 



r a («) = e- n °e -=i" +1 , (78) 

where the separation constant now takes the form f2 in which we absorb the 
Schwarzschild mass mo 



a = m n. (79) 

Consider, now, the case in which the background is not radiating. It is clear 
either from the solution above or from the original differential equation that for 
such a case the solution reduces to 



Y 2 (v s ) = e- n ° v °, (80) 

where the subscript s indicates quantities associated with the Schwarzschild 
geometry. One notices that in such a static background the quantity Y 2 (v s ) 
above constitutes the only time dependent part of the held. It follows then 
that to be consistent with the known [9 solutions we should require 



fl s = iuj, (81) 

where u is the frequency of the gravitational waves. This suggests that in the 
case of the radiating background we should expect the parameter f2 to be a 
complex function of A and u> such that 

Umfi(A,a;)=w. (82) 

The integration of the differential equation for Y^ 2 (v) follows the same trend 
and we find that 

„ -r y ^-a"d" +1 
y_ 2 (v) = e- rv e , (83) 

where 
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7 = m r, 



(84) 



r being a complex function of A and u) such that 

lim r = iu (85) 

A^O 

As would be expected from the theory of deferential equations the separation 
constants f2 and T can not have unique values. The individual solutions we 
obtain will therefore be representatives of classes of solutions. The range of these 
solutions is described in terms of the frequency spectrum of the gravitational 
field which, in our classical treatment, takes on continuous values. It will, later, 
be shown that by using certain conditions on the solutions the functional form 
of these separation constants can be more rigidly fixed. 

6.2 The functions X p (£) 

Following the integrations of the first order differential equations for Y p we are 
now left with the two equations for X 2 and X_ 2 to solve. These are respectively, 

e + 4Amo) ^ - C [e 2 - 4£ - 4m (4A - 0)] ^ + (e, £ - 20m O) X 2 (£) = 



(86) 



and 



e (e-Z + 4m A) ^ - { [5^ - 4£ + 4m r] ^ 



+ (C 2 + 2e t £ + 4m r) X_ 2 - 0. (87) 

It is clear that at £ = (or r = oo) both the equations above have regular sin- 
gularities [[l6]. This encourages us to seek for analytic solutions. Such solutions 
at £ = (r = oo) should be useful in discussing the asymptotic fall-offs of the 
fields and the question of energy flux. 

6.2.1 The peeling behavior 

Our initial goal is to develop asymptotic solutions for the functions X 2 (£) and 
X-2 (£)■ Consider a zero rest mass spin-s field tp p in a helicity state p. According 
to the peeling theorem by Roger Penrose E|, the quantities A s+p+1 ' > 'ip p and 
r (s-p+i)-0^ have a limit at null- infinity. In the case of gravitational fields we 
expect the outgoing components of the solutions to fall off as 

^Cp=±2) ~ r (s+ P +i) = r (2±2+i) ' ( 88 ) 
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while the ingoing solutions should fall off as 



VV=±2) ~ r(fl _ p +l) - r (2 T 2+l) ■ ( 89 ) 

It is necessary, therefore, that the solutions to our differential equations display 
the above asymptotic behavior. This, indeed, will be one of the tests for their 
validity. 

6.2.2 The Indicial Equations. 

It has been pointed out that at £ = we have a regular singularity in both 



equations 86 and |87] . Therefore it seems natural to attempt developing solutions 
about this point. Such solutions will be valid at far distances from the black 
hole. This class of solutions at such distances is useful if one is to engage, as we 
shall later, in a meaningful discussion of the gravitational energy flux. 

Let us assume that equation ^6] admits, as a solution, a series expansion 
about £ = of the form 



X 2 (?) = J2a n C +k , (90) 

n=0 

where k is some value to be determined. Using equation ^ in equation p6| gives 



Y^a n (n + k)(n + k-2) C +k " 2 

n=0 



a n [{n + k) (n + k - 5) - e{\ f 



+ 4a„ {(n + k)[(n + k — 5) Arn - m O] + 5rn ft} C +k = (91) 

n=0 

For n = 0, ao ^ we get the indicial equation 

k[{k-5)X-n]+5Q = 0. (92) 
which has two distinct roots, 



5,-). (93) 
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Similarly, for X^2 we can assume a solution of the form 



x- 2 (o = j2b n e +s , 



n=0 



(94) 



where, again s is some value to be determined. Substituting equation ^4] into 
equation 87, we find that 



^6„_ 2 [(n + s - 2) (n + s - 8) + 1] £ 



n+s 



-X>„-i [(n + « - 1) (n + « - 6) + 2e ; ] £" +s 



+4m ^fe„ [A (n + s) (n + s - 1) - T (n + s) + T] = (95) 

n=0 

From the equation above and the condition that bo =/= we obtain, for n = 0, 
the indicial equation 



r\ r 

A S+ A=° 



(96) 



whose roots arc 



1, 



A 



(97) 



Equations |93| and |97] (and the fact that Q and T are complex quantities) indicate 
that we can expect two linearly independent solutions for each of the fields. 

Asymptotic conditions The solutions to the indicial equations given in 
equations ^ and |97] fix for us the leading terms for the functions X2 (£) and 
X-2 (£), respectively. Thus from equation |93| we see that 



(98) 



and from equation Wi 



X_ 2 (0 ~ £ or ^ 



(99) 
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Both equations and show that the first solutions are consistent with the 
peeling theorem and can, in fact be recognized as outgoing fields (recall £ = 

2m(u) 2mo(l — Xv)\ 

r r ' ' 

On the other hand the second solutions are scaled by the quantities ft and 
r, respectively. These are the same arbitrary separation constants which, in 
the last section, we showed to be complex. Since, physically our solutions 
represent gravitational fields these constants must now be chosen to conform 
with the known boundary values for such ingoing waves. Consequently in order 
to satisfy the peeling theorem, it is clear that we must have Rett ~ A, so that 
X 2 ~ - and ReT ~ 5A, so that X^ 2 ~ -r- Moreover, the imaginary parts of 
these quantities must reduce to the limiting cases, lim D, — iui and lim Y = iu 

A^O A^O 

as was shown to be the case. These two conditions dictate that we set 

n = x + iuj, (ioo) 

and 

T = 5X + iuj, (101) 
The roots to the indicial equations |9^ and [)?], respectively now become 



5, 1 



(102) 



and 



1, 5 



(103) 



So that as £ — * 0, 



X 2 ^)^eorX 2 (0^^ u , (104) 

and 

X-2 (0 -» £ or X_ 2 (0 -» £ 5 £>. (105) 

The full functions X 2 (£) are readily obtained by writing down recurrence rela- 
tions using equations |9l| and The general solutions are, in each case, found 
to be linear combinations of the outgoing component X p OMt ^ and the ingoing 
component X p n \ 
Thus 



X p (0 = A p X^ + B p X^ . (106) 

Here the A p and B p are arbitrary constants of integration. These solutions are 
also found to converge. 
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7 The asymptotic solutions and physical infor- 
mation. 

7.1 Significance 

The principal aim of our study is to understand how gravitational waves are 
scattered by a background radiating spacctime. In particular, we are interested 
in the measurable physical results of this process, such as the energy flux and the 
manner in which the waves are reflected and absorbed by a radiating black hole. 
To this end we have, in the preceding discussions, developed field equations that 
describe the effects of these waves on the background spacetime. The physical 
quantities that we seek should, in principle, be calculated from the solutions of 
these equations. As can easily be shown, however, the series solutions obtainable 
are a result, in each case, of a three term recursion relation and so contain 
various coefficients that are not easy to relate. This feature of our solutions 
would seem to make inconvenient, their use in calculating a number of other 
physical quantities. It turns out, though, that for the features of our interest it 
is sufficient to consider the form of the solutions at certain special points. For 
example Chandrasekhar shows that a knowledge of the incident and reflected 
wave amplitudes can be deduced from the form of the solution at null-infinity. 
Moreover, one can also engage in a meaningful discussion pertaining to energy 
flux at these points. This means that we need only consider the leading terms 
in the solutions. 



7.2 The Source Terms 

In creating the function \&4 we have, so far, only considered solutions for the 
homogeneous part of the original differential equation (see equations [56| and 
[57]) . However, the full equation for ^4 is inhomogeneous so that the complete 
solution should contain a contribution due to the sources. We recall that the 
source term is scaled by the luminosity L = = Am which obviously 

vanishes as the background radiation is switched off, A — > 0. Since, in the first 
place, A < 1 => < 1, we shall presently assume that at large r values 

the source terms do not contribute significantly to the solution. Consequently 
we shall consider the asymptotic solutions from the homogeneous equation to 
be a sufficient representation of the general asymptotic solutions. With this we 
now write down the asymptotic form of the entire solutions. 

7.3 The Solutions 

It was shown, in equation [78|, that for vto, 



Y 2 (u) = e-xMi-A.) = e -nu e ^A^" v " + \ ( 107 ) 
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For small A values, the expression for the logarithmic expansion can be written 
to first order in A. This gives 



Y2(u) 



ln(l-Au) ^ g-^ 1 -' 



(108) 



And using the conditions spelled in equation 10C to satisfy the peeling property 
for the ingoing field, we find that 



Y 2 (u) ~ e - Av -™ v (109) 
Similarly, going through the same treatment for Y_2 (u) and applying equation 



101, we find that 



Yn (u) 



-5Xv — iujv 



(110) 



The above results along with the angular solutions pY™ (9, ip) = P S™ (9) e lmtp 



of equation |60j can now be joined to the functions of equations 104 and 105 
(recall £ = m ^°' — 2m "(i~ A ") ) t gj ve the following asymptotic solutions for 
the functions and ^4- 



outgoing 



2Yr(oM 



[2m(u)]° 



ingoing 



2m{u) 



ym N")] S e -B(^)« e ip e -ato« 

(111) 



*4 ~ _ 2 Y l m {9,cp)^^e- 5 (^) v e-^ v 

Here, p(r) = j In (^ a ) and where for physical reasons we find it useful to 
express the solutions in terms of the luminosity L as given by L = — = Atoo- 

8 Conclusion 

We have obtained analytic solutions to the problem of gravitational fields propa- 
gating in a radiating spacetime. These solutions satisfy all the known conditions 
for the propagation of spin-s, zero rest mass fields. From their asymptotic form 
in equations 111 it is seen that the solutions are completely consistent with the 



peeling theorem of Penrose and fall off in the manner predicted by this theorem. 
Moreover, one observes, further, that as the background radiation is switched 
off (i.e. in the limit A — > 0), the theory recovers the known solutions (see for 
example ||) for the perturbed static geometry of Schwarzschild. We consider 
the passing of these two tests a validation of our analysis. 

One new significant feature this analysis brings to surface is that the so- 
lutions (see equation 111) are scaled by factors of the form e ^ +P+1 H "»o) u ) 
where s = 2 and p — ±2. But is positive definite. Consequently, these 
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factors indicate that when gravitational fields propagate in a radiating space- 
time they suffer an attenuation, and this attenuation can be quantitatively 
described. The attenuation weight seems in turn to be directly related to the 
spin weight of the perturbed fields. It is also scaled by the luminosity L of the 
background. Further, as one notices from the solutions, the attenuation per- 
sists independent of whether the fields are ingoing or outgoing. It is of course 
fair to ask whether this character of our solutions is not, in the first place, a 
reflection of the mass function that we chose. Recalling that general radiative 
mass function to (v) is a monotonic decreasing function in v an expansion of the 
m (v) = too — Lrrif) — ^J^too — .... about v = indicates that the first order term 
in the luminosity would seem to make the significant contribution. This seems 
to suggest that the attenuations manifested in our solutions are independent of 
the manner in which the blackhole radiates and may persist for any mass func- 
tion chosen. As far as we know this seems to be a new feature in the literature 
of this branch of general relativity; one that may, indeed, have some interesting 
astrophysical implications. 

A persistent attenuation of this sort would seem to suggest the possibility 
that energy is being dumped into the host spacetime. Over large time scales, 
this could have significant implications on the evolution of such a radiating 
system. This question can, however, only be resolved by a rigorous calculation 
of the energy flux. For such a calculation and an extension of this discussion 
see |i~8| We intend to follow up this issue in future discussions. 
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